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Topological Josephson junctions carry 47r-periodic bound states. A finite bias applied to the 
junction limits the lifetime of the bound state by dynamically coupling it to the continuum. That 
lifetime may be shorter or longer than the phase adjustment time for a given bound state, depending 
on the resistance of the circuit "seen" by the junction. We show that, in the former case, the presence 
of a 47r-periodic bound state manifests itself in a peak in the current noise spectrum at half the 
Josephson frequency, when a constant bias is applied to the junction. In the opposite case of fast 
phase adjustment, the 4-7r-periodicity is evidenced by an even-odd effect in Shapiro steps. We specify 
the conditions necessary for observing the manifestations of 47r-periodicity in the noise spectrum 
and Shapiro step measurements. 

PACS numbers: 71.10.Pm, 74.45. +c, 05.40.Ca, 03.67.Lx 



Topological Josephson junctions have attracted much 
interest lately as a means of probing the zero-energy Ma- 
jorana fermion states that exist at the surface of topolog- 
ical superconductors. Such topological superconductors 
may be realized via the proximity effect by combining 
conventional superconductors with two-dimensional (2D) 
topological insulators [l| or with nanowires in the pres- 
ence of both strong spin-orbit coupling and a magnetic 
field Recently several experiments have reported 

evidence of zero-energy states in nanowirc-based systems 
[3-0] ■ To confirm their Majorana nature, additional ex- 
perimental signatures are desirable. 

In a topological Josephson junction, the Majorana 
bound states localized on either side of the junction hy- 
bridize and form an Andreev bound state whose energy 
ca(v) is 47r-periodic in the phase difference <p between the 
two superconductors. Depending whether the state is oc- 
cupied or empty, the energy of the junction is ±£a(}p)/2. 
In the presence of parity-changing processes, the occu- 
pation of the state may change. Thus, the equilibrium 
Josephson current displays the usual 27r-periodicity as 
the system follows the ground state. By contrast, if 
upon phase variation the system follows one branch of 
the spectrum, then 47r-pcriodicity should appear indeed. 
As a result, under dc bias voltage Vd c , such a system has 
beenpredicted to manifest a fractional ac Josephson ef- 
fect [3,0, Hi at frequency uij/2 — eVd c /ft, that is at half 
of the "usual" Josephson frequency. By the same token, 
in the presence of an additional ac bias with frequency 
CI, one would expect an even-odd effect: namely only 
the even Shapiro steps at eVdc = kCl (k g Z) should be 
visible in the current- voltage characteristics [8r|l0| . How- 
ever, the application of a bias voltage inevitably couples 
the bound state to the continuum, thus causing its occu- 
pation to switch. The corresponding switching rate de- 
termines the lifetime of the bound state, r s . In addition 



to these intrinsic processes, the evolution of the phase 
difference across the junction depends on the properties 
of the circuit connecting the Josephson junction to the 
voltage source. Any non-zero resistance 1Z of the connec- 
tion allows for an adjustment of the phase difference over 
some characteristic time [ll], E3 T R ^ T^ 1 - 

In this work, we evaluate the lifetime of Majorana 
bound states, r s , limited by their dynamic coupling to the 
continuum. This mechanism gains importance in nearly- 
ballistic junctions and leads to a strong dependence of 
r s on the applied voltage. We show that the transport 
properties of the junction crucially depend on two char- 
acteristic time scales, r s and tr. If r s 3> Tr, Majorana 
states lead to an even-odd effect in the height of Shapiro 
steps, in agreement with Rcfs. By contrast, if 

t s <C tr, all Shapiro steps are suppressed. However, sig- 
natures of the 47r-periodicity are still visible in the noise 
spectrum which, under dc voltage bias, displays peaks at 
lu = ±luj/2, as was seen in numerical simulations [l3| . 
Here we develop an analytical theory for the noise spec- 
trum and find the dependence of the peak widths on t s . 
While noise measurements in the GHz-range are not easy 
to realize, the low-frequency noise is more accessible. The 
down-conversion of the noise peak to w = ±(lo,j/2 — CI) 
may be achieved by adding a small ac bias of frequency 
Cl. 

To examine the non-adiabatic transitions between the 
Majorana state and quasiparticlcs continuum, we con- 
sider the helical edge state of a 2D topological insulator 
in which superconductivity has been induced by two su- 
perconducting contacts in order to create a topological 
Josephson junction of length L. The system is described 
by the Hamiltonian 

H = vpa z r z - eU(x,t)r z + M(x)a x + A(x)e t,i,ix ' t)T 'T x . 

(1) 

Here, v is the Fermi velocity, p is the momentum opera- 
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tor, U(x,t) = V(t)[9(-x)-9(x-L)]/2 is the electric po- 
tential, M(x) = M6(x)9(L — x) is a transverse magnetic 
field within the junction, A(x) = A[6>(— x)+6*(a- — L)] and 
cj)(x,t) = ip(t)[9(-x) - 9{x - L)]/2, with y>(i) = 2eV(t), 
are the amplitude and phase of the superconducting or- 
der parameter in the left and right leads, and er,;, tj 
= x,y,z) are Pauli matrices acting in the spin and 
particle/hole spaces, respectively. All energies are mea- 
sured from the chemical potential. 

We concentrate on the case of a short junction, L -C £, 
where £ = v/A is the superconducting coherence length. 
In equilibrium (V — 0), such a junction hosts a single 
Andrecv bound state with energy 

e A (if) = VDAcos(<p/2), (2) 

where D is the transmission probability of the junction 
which depends on its length and on the magnitude of the 
transverse field. Thus, the minimal gap S between the 
bound state and the continuum at <p = 2mt (n £ Z) is 
given as S = A(l — \/~D). In the following we consider a 
highly transmitting junction where 5 f=a AR/2 and R = 
1 — D ~ (ML/v) 2 is the reflection probability. 

Out of equilibrium, non-adiabatic transitions between 
the Andreev bound state and the continuum are induced. 
These transitions change the occupation of the bound 
state and thus lead to switching between the two current 
branches, I(p) = ±Jj sin(<p/2), where Ij = e\fT) A/2. 
At dc bias eVdc *C A, violation of the adiabaticity oc- 
curs in narrow intervals \cp — 2ixn\ <C 7r of the time- 
varying phase if = 2eVd c t- To find the correspond- 
ing probability of a non-adiabatic transition between 
the localized Majorana state and continuum, we con- 
centrate on the case n = 0, corresponding to the time 
interval \t\ <C 7r/(eVd c ). Using a gauge transforma- 
tion H -> H = WHU - iWU with U = exp[i0T 2 /2] 
and taking the limit L — >• (keeping i? fixed), we ob- 
tain for the said interval of ip the simplified Hamiltonian 
H = vpa z r z + At x + v[(tp/2)a z + VRij x ]5(x). Further- 
more, at eVdc *C A, only states close to the continuum 
edge, v\p\ -C A, are relevant. Thus, after diagonalizing 
the bulk Hamiltonian, we can restrict ourselves to a 2 x 2 
subspacc of the initial spin and particlc/hole space, 

H = A + ^- + v^a z + VRa^j6(x). (3) 

Eq. (|3]) describes a spin-degenerate continuum with 
quadratic dispersion, in the presence of a spin-dependent 
local potential. The first term in this potential accounts 
for the phase shift across the barrier in a gauge with zero 
electric potential in the leads and a vector potential lo- 
calized at the barrier. The second term describes the 
magnetic barrier. 

For a fixed phase, the Hamiltonian ([3]) accommo- 
dates a single bound state with energy €a(<p) — 
A (l - ip 2 /8 - R/2), in agreement with Eq. © at 



R, ip 2 <C 1. A particle occupying this bound state at 
time t — oo (within the simplified model) has a proba- 
bility s to escape to the continuum as the phase increases. 
The problem is, thus, a generalization to a two-band 
model of the transition from a discrete state to a contin- 
uum, considered by Demkov and Osherov [lil ]. Dimen- 
sional analysis shows that the transition (or switching) 
probability is determined by the adiabaticity parameter 
A = i? 3 / 2 A/(eVdc). Below we find this probability in two 
limiting cases of the parameter A. 

Let us start with the anti-adiabatic regime, A C 1. 
At A = the spin bands in Eq. ([3]) are decoupled. At 
times t < 0, the bound state belongs to the spin-up band 
whereas, at times t > 0, the bound state belongs to the 
spin-down band. The spin-up bound state is described by 
a wavefunction \ip^(t)) whose projection on the position 
of the local potential is 

(x = 0|VtW) = -4= [ du e-' e *(~ 2 -) 3/2 /3| t ). (4) 
V 2ttI Jc 

The wave function \tp±(t)) of the spin-down bound state is 
related to |^f(i)) by time reversal. Here, the characteris- 
tic length and time scales are given by I = w/[A 2 eVd c ] 1 ^ 3 
and t = l/[A(eVdc) 2 ] 1 ^ 3 , respectively. Furthermore, C is 
a contour in the complex w-plane [14J that starts and 
ends at infinity, with arguments 7r < 8 < 57r/3 and 
< 9 < 7r/3, respectively, and avoids the branch cut 
along the positive real axis. As the spin bands are de- 
coupled, a particle occupying the (spin-up) bound state 
at t — —oo, has probability 1 — s = to occupy the 
(spin-down) bound state as t — > oo. 

A finite A couples the two bands and, thus, enables spin 
flips. The switching probability s can be obtained from 
the overlap c±{t) — {ip±(t)\ij}(t)) of the exact wavefunc- 
tion, \ip(t)), where \if>(— oo)) = \ip^(—oo)) , with the wave- 
function of the spin-down bound state, \ipi(t))-i through 
s = 1 — |cj.(oo)| 2 . At A <C 1, cj. can be computed pertur- 
batively using 

c^t) a i(^(t)\vVRa x 6(x)\^ t (t)). (5) 

Solving the differential equation ([5]) to obtain 04.(00) and 
computing s, we find s « 1 — 1.05 A 2 / 3 . The time scale 
over which the transition happens is r t ~ r. 

In the quasi-adiabatic regime, A 3> 1, it is convenient 
to expand the exact wavefunction in the adiabatic basis 
of Eq. ©, 

|V(t)) =c A {tM A (t)) +X) C ^(*)IV(*))- ( 6 ) 

pa 

Here \if>A(t)) and \i(> p< r{t)) are the adiabatic wavefunctions 
for the bound state and the doubly degenerate states of 
the continuum, respectively. [Note that \iPa(To°)) = 
1^4(^00)).] The switching probability s is related to 
the amplitudes c pcr (oo) of the continuum states in Eq. © 
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FIG. 1: Switching probability s as a function of the adiabatic- 
ity parameter A (dots). The lines correspond to the asymp- 
totic expressions for s, see text. 
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FIG. 2: Schematic view of the switching processes due to the 
coupling with the continuum: (a) occupation n of the bound 
state, (b) energy e of the system, and (c) Josephson current 
7 as a function of time t under dc bias voltage Vdc. 



through s = X)pcr l c po-(o°)| 2 j using the initial conditions 
Ca{— oo) = 1 and c pa (— oo) =0. At A 1, using 

e p( T — eA[t) 

(7) 

the amplitudes c pa (oo) are expressed through integrals 
that may be evaluated by a saddle point method. We 
obtain s ~ 0.93 A -5 / 4 e -2A / 3 . Furthermore, we can iden- 
tify the time scale over which the transition happens, 
r t ~ VR/(eV dc ). 

At arbitrary A, the switching probability can be 
obtained numerically by discrctizing Eq. (J3j) on a 
tight-binding lattice and solving the corresponding 
Schrodinger equation numerically. The result, together 
with the asymptotes obtained above, is shown in Fig. [TJ 

Using the fact that the transition time rt is much 
shorter than the Josephson oscillation period, r t <C 
7r/ (eVdc), we may now write an effective discrete Markov 
model for the bound state dynamics, cf. Fig. [U Using 
a discrete time evolution we assume that if the state is 
filled, at phase ip2n = 4n7r, there is a probability s of the 
particle to escape from the bound state to the continuum, 
whereas if the state is empty, at time tp2n+i = (4n + 2)7r, 
there is a probability s of a particle from the continuum 
filling the bound state. Thus, 



Pra 
Qln 



1 S \ P 2n 

l-sj \Q 2n 

l-s 0\ /P 2 „_! 

s 1/ \Q2n-1 



(8b) 



Here P n is the probability for the state to be occupied, 
and Q n = 1 — Pn is the probability for the state to be 
empty at phases ip n < ip(t) < (p n +i, corresponding to 
n = lnt[ip(t)/(2iv)]. Solving these equations iteratively, 
we obtain 

P2(n + k) - ^ + (l-s) 2fc (P 2)l -P-), (9a) 

P 2 (, i+fc ) + l = P^ + l + a-s) 2k+1 (P2n~P 2 00 n ). (9b) 



The long-time probabilities at k ^> — l/ln(l — s), cor- 
responding to t S> t s — — 27r/[eVd c m(l — s)}, are 47r- 
periodic and independent of the initial state: P^ = 
[l-(-l)V(2-s)]/2. 

In order to determine the transport properties of the 
junction, the switching time r s has to be compared to 
other characteristic time scales of the system. In par- 
ticular, if the junction is embedded into a circuit with 
a resistance 1Z in series, the phase difference across the 
junction may adjust over a typical time scale tr oc Tlr 1 . 

If t s 3> tr, switching may be neglected (i.e., the 
current may be obtained using the initial occupations 
Pq/Qo). Computing the dc current in the presence of an 
applied voltage V(t) = Vd c + U ac cos(f2£) with V ac <§; Vd c 
then yields the even-odd effect discussed in Rcfs. 
Namely, taking into account a finite resistance 1Z, the 
average current reads 




(— V 



( TUk 

\sv k 



(10) 

where I k = Ij\Jk{(x)\ is the height of the Shapiro step 
at eVdc = kVt, and 8Vk = Va c — kCl/e. Here J k are the 
Bessel functions and a = eV ac /f2. The characteristic 
time scale may be identified as [H, IH 4 k) = l/(eKI k ) 
at eV^c ~ kfl. 

A small resistance satisfying the relation TZIj 17 is 
advantageous for the resolution of Shapiro steps. Fur- 
thermore, long switching times require s <C 1. In 
the small- s regime, the time scale t s decreases as 
cxp[2P 3 / 2 A/(3fcfi)] with increasing k. On the other 
hand, at small ac perturbation, a <C 1, the time scale 
rj^ increases exponentially with k. The crossover from 
T s 3> t~r to the opposite limit, t s <C tr, upon increasing 
Vdc may occur without violation of the condition s«l. 
This restricts the number of observable Shapiro steps in 
the current-voltage characteristics, as we show now. At 
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t s <C Tr, we may use the long-time probabilities P°° /Q°° 
to compute the current and take the limit tr — > oo. At 
long times, the average current is 27r-periodic, 



(7(i)) = Ij sin 



V(*) 



(Qln1#(t)/(27r)] 



300 \ 
Int[ V (i)/(2 7 r)]/ 

(11) 



More importantly, Eq. (|lip shows that the current is pro- 
portional to the switching probability s <C 1. This result 
is valid in the presence of microwave irradiation as long as 
V ac <C Vdtn O/e , su ch that the phase velocity (p remains 
always positive [L7f . Thus Shapiro steps are strongly sup- 
pressed, (Ik) oc s. 

In order to reveal signatures of the 47r-periodicity in 
the regime r s -C tr, we now turn to the current noise 
spectrum, 



S(u) = 2 dTcoa(urr)(8I(t)6I(t + T)), (12) 



where SI = I — (I) and the bar denotes time-averaging. 
It may be obtained from the correlator (I((pi)I(ip2)) = 
1} 8111(^/2) sm(<p 2 /2)[Q™x n2 (P ni =0) - P™x n2 (P ni = 
1)] at ip 1 < <f2, where = Int[y>j/(27r)]. Using the 
conditional probabilities obtained from Eqs. ([9|), we find 



(61^)61^2)) 



4/3(1 



sin — sin— (1 - s)™ 2 "™ 1 . 
s) 2 2 2 X ' 

(13) 

At dc bias only, the noise spectral density evaluates to 



(2 



4slj 



(eV dc f 



4 cos' 



2eV dc 



n(2-s) [ W MeT/ d c) 2 ] 2 4cos 



If s <C 1, it has sharp peaks at lu -- 
the "usual" Josephson frequency: 



2eV dc 1 

±eV dc , i-e. 



(14) 



at half of 



S(u) 



I 2 j 



seV dc /ir 



2 {lu t eV Ac f + (seV dc /iry 



(15) 



at \lo =f eVdd <C eV dc . In particular, the peak width 
is 2seV dc /ir. The position of the peak reveals the 47r- 
periodicity of the Andrcev bound state whereas the in- 
verse width characterizes its lifetime r s oc s _1 . The peak 
in the noise is due to the transient 47r-periodic behavior 
18| of the current at times smaller than the lifetime of 
the bound state. 

Under microwave irradiation, the peak may be shifted 
to smaller frequencies. In particular, in the limit U ac <C 
V dc , fi/e, we find 



seV dc /ir 



s{u) - Y Jk{a) \^leV dc -kn)}2 + ( se u dc /7r)2 



(16) 



at |w=F (eVd c — <C eV dc . As above, the peak width is 
set by the lifetime of the bound state which, thus, may 



be probed by noise measurements. Eq. (fT6|) holds for 
frequencies to not too close to zero. In the limit to 
0, additional features related to the Shapiro steps may 
appear [l6| . 

While we considered the helical edge states of 2D 
topological insulators, the model is also applicable to 
nanowires [H, 0| with strong spin-orbit coupling and a 
Zeeman energy much larger than A. Note that, in addi- 
tion to the non-adiabatic processes that we considered, 
non-adiabatic processes in the vicinity of ip = (2n + 1)tt 
become important if the zero-energy crossing is split 
due to the presence of additional Majorana modes at 
the ends of the wire 18 -2^. In particular, in order 
to see signatures associated with the 4-7r-periodicity, the 
probability of Landau- Zener tunneling across the gap at 
ip = (2n + 1)tt would have to be large while the switch- 
ing probability due to the coupling with the continuum, 
discussed in this work, remains small. 

To summarize, we studied the characteristics of the 
voltage-biased topological Josephson junctions. The dy- 
namic coupling between the Andrcev bound state and the 
continuum yields a finite lifetime for the bound state. If 
that lifetime is longer than the characteristic time scale 
over which the phase in the circuit evolves to gener- 
ate Shapiro steps, these steps display an even-odd ef- 
fect. However, if the lifetime becomes comparable to that 
timescale, signatures in the average current are obscured 
while the current noise remains sensitive to the under- 
lying 47r-periodicity of the bound state, which leads to 
characteristic peaks in the noise spectrum. The peaks in 
the finite-frequency noise, at dc bias, and low-frequency 
noise, with additional ac bias, should be seen easily at 
voltages eVdc < i? 3 ^ 2 A, where R is the reflection proba- 
bility. The peak width is a measure for the rate of parity- 
changing processes. 

In the final stages of preparing the manuscript, we be- 



came aware of Ref. [21[ considering related effects in 



nanowire-based topological Josephson junctions. 
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SUPPLEMENTAL MATERIAL 



The RSJ-model for the topological Josephson junction 

The voltage-biased RSJ-model [l| is described by the equation 



V(t) = KI s (t) + -<p{t), 



where V(t) = V dc + V ac cos (fit) and I s (t) = Ij sin [ip{t)/2]. 
We rewrite 



yielding 



tp(t) = 2knt 



eV dc - kfl = elZIj J - 



2eV ai 



neV a < 

n 



sm{nt) + X (t), 



(fc-m) nt + - x (t) 



(17) 



(18) 



Keeping only the slowly varying contribution m = k 7 we find 



eVdc — kfl ~ elZIjJ- 



(19) 



for |eVdc — kfl\ <C eVd c - For A; = 0, Eq. (|T9|) describes the supercurrent branch at low bias. The Shapiro steps 
at eVdc = kQ arc replicas of the supercurrent branch Q with a reduced maximal current 1^ = Ij\Jf.(a)\, where 
a = eV ac /n. 

From Eq. (fT9|) . we extract the characteristic time scale rjj = l/[elllj\ Jk(ce)\] = l/(elZIk)- For |(eVdc — kfl)r^ \ < 

1, we find the constant solution sinx = (eVdc — fcfi)r^ , whereas, for |(eVd c — fcfl)rjj | > 1, integration of Eq. (fl9|) 
over one period T yields 



T 



2ttt 



ik) 



(20) 



[eV dc - 



(fe) 



The dc current I dc = Is is then given as 



idc - IjJ-k (a) sin— — - 



1 

eft 



eU dc - kQ - -x(t) }• , 



(21) 



yielding Eq. (10) in the main text. 



Long-time properties of the topological Josephson junction in the limit tr — > oo 

In the presence of switching such that t s <C tr, we can get analytic results for the situation of perfect voltage bias, 
11 = 0. Using the long-time probabilities, cf. Eqs (9) in the main text, we obtain 

</(<)> = ^ Ij Y~s ^ ^ J- k (2n-^\ C os[{2neV dc -m)t + n Vo ], (22) 

n — — oo fc— — oo 

where </?o is the initial phase. 

In particular, the dc current reads 

(-fee) = ^2^7 S i _4 n 2 £ ^-fc(2n^^cos(n^o)5 2ne ^ cjfcn , (23) 

n— — oo A;— — oo ^ ' 
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consisting of a "background" current independent on the initial phase, 



(I 



(bg) s 
dc / 



2 8 



(24) 



as well as delta functions at 2neVd c = kQ for n ^ 0, with an amplitude dependent on the initial phase. 
Similarly, we may also compute the finite-frequency noise, 



E 



J_; [(2m + 1 — ix)a] J-(i+k) [(2(m + n) + 1 — 2ix)a] 



e r 02neV dc ,kU 



4ir 2 (2 — s) 2 ^—f (m — ix)(m + 1 — ix)(m + n — ix)(m + n + 1 — ix) (2m + 1 — ia;)ey — =F <*> 



(25) 



Here the "background" noise independent on the initial phase reads 



4tt 2 (2 



E 

mj;± 



J_ ( ((2m + 1 - ix)a) 
(m — ix)(m + 1 — ix) 



1 



(2m + 1 — ix)eVd c — ttl^f uj 



(26) 



For s <C 1, the noise S^ bg '(u}) has peaks at w = ±(eVd c — Jfi). Namely, 

5 (hff) (w 



2 J / l ^WT^c-ZO^ + ^dc/TT) 2 ' 



(27) 



for |cj =F (eVdc — *C eVd c - As for the dc current, in addition, there are delta functions at 2neVd c = k£l for n ^ 0, 
with an amplitude dependent on the initial phase. 

The characteristic time scale for the current and noise to become independent of the initial phase is tr. Thus, to 
obtain the shape of Shapiro steps as well as of the additional features in the noise spectrum, it would be necessary 
to determine the occupation probabilities P(ip,t) and Q(ip,t) in the presence of a finite series resistance TZ. We leave 
this for further study. 
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